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Selection  of  Near  Neighbors  in  VDF 
Limits  Thermalization 


Merge 


E(n-\-  2) 
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Review  of 


Merge  to  Pair  — »  DOF  for  Conservation: 

o  (n+2):2  yields  Exact  Mass, 
Momentum,  and  Kinetic  Energy 
Conservation 

o  Applied  Spatially  also  Shown  to 
Conserve  Electrostatic  Energy 

©  Though  Energy  Conserving, 

Still  Thermalizes  VDF 

Selection  of  Near  Neighbors  in  VDF 
Limits  Thermalization 

Merge  via  Separate  Octree/Species 
Only  Change  for  Mixtures ! 


Distribution  A:  Approved  i 
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Octree  Velocity  Bins 


Efficient  Neighbor  Selection 


From  RGD30:  Mach  2 


ID  Normal  Argon  Shock  Test 
o  Simple  Verification  vs.  DS IV 

o  Initial  Conditions: 

To  =  293K,  n0  =  lE22/m3,  v0  =  637.4(m/s) 
o  Initial  Jump  to  Post-Shock  at  1cm 

©  VHS  Collisions: 

T ref =21 3K,  dref=4.HA,  CJy//5=0.81 
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o  Initial  Conditions: 
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©  Time  Average: 
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From  RGD30:  Mach  2 


ID  Normal  Argon  Shock  Test 

o  Simple  Verification  vs.  DS IV 

o  Initial  Conditions: 

T0  =  293K,  n0  =  lE22/m3,  v0  =  637.4(m/s) 

o  Initial  Jump  to  Post-Shock  at  1cm 

o  VHS  Collisions: 

Tref=213K,  dref=4.llA,  ccv//,s'=0.8 1 

©  Time  Average: 
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©  Error  Controlled:  err  oc  v/ N/cell 
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2D  Argon  Shock  Test 

o  Initial  Conditions  like  M=2  Except: 
vo  =  2550m/s 

©  Specular:  x=5  —  5.04mm  with  y=±2mm 

©  Half  Domain  Modeled: 

80/xm  x  80/xm  Cells 


gon  Bow  Shock 


TURF:  n  -  Standard  DSMC 


X  (mm) 

TURF:  n  -  SWPM+Octree 


11/21 


From  RGD30:  Mach  8 


2D  Argon  Shock  Test 

o  Initial  Conditions  like  M=2  Except: 
vo  =  2550m/s 

o  Specular:  x=5  —  5.04mm  with  y=±2mm 
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From  RGD30:  Mach  8  Argon  Bow  Shock  o 


2D  Argon  Shock  Test 

o  Initial  Conditions  like  M=2  Except: 
vo  =  2550m/s 

o  Specular:  x=5  —  5.04mm  with  y=±2mm 
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©  Time  Average: 
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o  SWPM  Similar  to  Standard  DSMC 
o  Despite  Different  Np /  Cell 
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TURF:  n  -  SWPM+Octree 
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X  (mm) 
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He:Xe  Mixture  Shock 


ID  Normal  He:Xe  Shock  Test 

o  Mach  3.89  with  He:Xe  of  97:3 
(i.e.  Bird ’94  Fig  12.35) 


Robert  Martin  (AFRL/RQRS) 


DlSTRIBU] 


Converged  lOOx  Baseline 


:  □ 


& 


'O  A  O' 


He:Xe  Mixture  Shock 


ID  Normal  He:Xe  Shock  Test 

o  Mach  3.89  with  He:Xe  of  97:3 
(i.e.  Bird ’94  Fig  12.35) 

©  Highlights  Species  Separation 
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He:Xe  Mixture  Shock 


ID  Normal  He:Xe  Shock  Test 

©  Mach  3.89  with  He:Xe  of  97:3 
(i.e.  Bird ’94  Fig  12.35) 

©  Highlights  Species  Separation 
a  Separation  Peak  at  pHe  ~  pxe 
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He:Xe  Mixture  Shock 


ID  Normal  He:Xe  Shock  Test 

©  Mach  3.89  with  He:Xe  of  97:3 
(i.e.  Bird ’94  Fig  12.35) 

©  Highlights  Species  Separation 
©  Separation  Peak  at  pHe  ~  pxe 

©  DSMC  needs  33x  He:Xe  Macroparticles 
(From  mXe/mHe ) 

©  Reduced  Particle  Count  Introduces  Error 
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He:Xe  Mixture  Shock 


ID  Normal  He:Xe  Shock  Test 

©  Mach  3.89  with  He:Xe  of  97:3 
(i.e.  Bird ’94  Fig  12.35) 

©  Highlights  Species  Separation 
©  Separation  Peak  at  pHe  ~  pxe 

©  DSMC  needs  33x  He:Xe  Macroparticles 
(From  mXe/mHe ) 

o  Reduced  Particle  Count  Introduces  Error 
©  Error  in  Time  Average  at  Baseline 


Robert  Martin  (AFRL/RQRS) 


DlSTRIBU] 


Comparison  Baseline  vs.  lOOx 


□  <i  g 


release;  Distribution  Unlimited;  PA  #17517 


^  'O  A  O' 


12/21 


He:Xe  Mixture  Shock 


ID  Normal  He:Xe  Shock  Test 

©  Mach  3.89  with  He:Xe  of  97:3 
(i.e.  Bird ’94  Fig  12.35) 

©  Highlights  Species  Separation 
©  Separation  Peak  at  pHe  ~  pxe 

©  DSMC  needs  33x  He:Xe  Macroparticles 
(From  mXe/mHe ) 

©  Reduced  Particle  Count  Introduces  Error 
©  Error  in  Time  Average  at  Baseline 
©  Instantaneous:  Dramatic  Noise  in  Xe 
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He:Xe  Mixture  Shock 


ID  Normal  He:Xe  Shock  Test 

©  Mach  3.89  with  He:Xe  of  97:3 
(i.e.  Bird ’94  Fig  12.35) 

©  Highlights  Species  Separation 
©  Separation  Peak  at  pHe  ~  pxe 

©  DSMC  needs  33x  He:Xe  Macroparticles 
(From  mXe/mHe ) 

o  Reduced  Particle  Count  Introduces  Error 
©  Error  in  Time  Average  at  Baseline 
©  Instantaneous:  Dramatic  Noise  in  Xe 

Noise  Reduction  via  Variable  Weights? 
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He:Xe  Mixture  Shock  -  Fractional  DSMC 


ID  He:Xe  Shock  with  SWPM+Octrees 
o  Xe  Noise  Controlled  by  1 : 1  Target 
®  He:Xe  Noise  Comparable 


lOOx  vs.  FDSMC  1:1 
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He:Xe  Mixture  Shock  -  Fractional  DSMC 


ID  He:Xe  Shock  with  SWPM+Octrees 
©  Xe  Noise  Controlled  by  1 : 1  Target 
©  He:Xe  Noise  Comparable 
©  Direct  Noise  Control  by  Target  Ratio 


lOOx  vs.  FDSMC  2:1 
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He:Xe  Mixture  Shock  -  Fractional  DSMC 


ID  He:Xe  Shock  with  SWPM+Octrees 
©  Xe  Noise  Controlled  by  1 : 1  Target 
©  He:Xe  Noise  Comparable 
©  Direct  Noise  Control  by  Target  Ratio 
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He:Xe  Mixture  Shock  -  Fractional  DSMC 


ID  He:Xe  Shock  with  SWPM+Octrees 
©  Xe  Noise  Controlled  by  1 : 1  Target 
©  He:Xe  Noise  Comparable 
©  Direct  Noise  Control  by  Target  Ratio 
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He:Xe  Mixture  Shock  -  Fractional  DSMC 


ID  He:Xe  Shock  with  SWPM+Octrees 
©  Xe  Noise  Controlled  by  1 : 1  Target 
©  He:Xe  Noise  Comparable 
©  Direct  Noise  Control  by  Target  Ratio 
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He:Xe  Mixture  Shock  -  Fractional  DSMC 


ID  He:Xe  Shock  with  SWPM+Octrees 
9  Xe  Noise  Controlled  by  1 : 1  Target 
9  He:Xe  Noise  Comparable 
9  Direct  Noise  Control  by  Target  Ratio 
9  Converged  Error  Finite!  (1:1  Target) 
o  Error  Source  Still  Unidentified... 
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He:Xe  Mixture  Shock  -  Fractional  DSMC 


ID  He:Xe  Shock  with  SWPM+Octrees 
9  Xe  Noise  Controlled  by  1 : 1  Target 
9  He:Xe  Noise  Comparable 
9  Direct  Noise  Control  by  Target  Ratio 
9  Converged  Error  Finite!  (1:1  Target) 
o  Error  Source  Still  Unidentified... 

9  Potentially  Sensitivity  to  He  Tails? 
-Merge  Impacts  Higher  Moments 
-TBD  Error  vs.  He-Noise  Level 
-Improvement  Merge  to  Preserve  Tails 

9  Adaptation  of  SWMP  Incorrect? 

-#  of  Collisions  Sampled  as  WDF  Varies? 
-Collision  Pair  Rule  Wrong,  wL  <  w7  ? 
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He:Xe  Mixture  Shock  -  Fractional  DSMC 


ID  He:Xe  Shock  with  SWPM+Octrees 
9  Xe  Noise  Controlled  by  1 : 1  Target 
9  He:Xe  Noise  Comparable 
9  Direct  Noise  Control  by  Target  Ratio 
9  Converged  Error  Finite!  (1:1  Target) 
o  Error  Source  Still  Unidentified... 

9  Potentially  Sensitivity  to  He  Tails? 
-Merge  Impacts  Higher  Moments 
-TBD  Error  vs.  He-Noise  Level 
-Improvement  Merge  to  Preserve  Tails 

9  Adaptation  of  SWMP  Incorrect? 

-#  of  Collisions  Sampled  as  WDF  Varies? 
-Collision  Pair  Rule  Wrong,  wL  <  w7  ? 

Error  Identification  in  Future  Work 
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Rothe’s  He:Ar  Free-Jet  Experiment 


o  Helium- Argon  Mixture 


Re=533 


o  Expanded  through  Nozzle  to  Vacuum 
o  e-Beam  Concentration  Measurements 

Jet  Flow  Layout 

r=96mm 


x=248mm 
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Prior  Results  from  DS2V  and  Continuum  ^ 


Case  Tested  in  PhD  &  RGD27 : 

Continuum  Boundary  Conditions 


dXAr  /dx=0 
dpv/dx=0 


dXAr/dx=0 

dpv/dx=0 


Monte  Carlo  Boundary  Conditions 

Surrogate  for  FRC  Injection 


Dynamic  Range  »  Shocks 


Solution  from  DS2V 


Solution  from  SPARTAN 
(Navier-Stokes  +  Diffusion  Velocity) 


Features  Found  via 


O  Wing  Increases  near  Nozzle  Edge 


Experimental 


DS2V 
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Features  Found  via  Simulations 


<  □  >  < 


<00,0 


Features  Found  via  Simulations 


O  Wing  Increases  near  Nozzle  Edge 
O  Lower  Radial  Boundary  Edge  Concentration 
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Features  Found  via  Simulations 


O  Wing  Increases  near  Nozzle  Edge 
9  Lower  Radial  Boundary  Edge  Concentration 
Q  Deeper  Jet  Edge  Concentration  Drop 


Challenges  for  Jet  in  TURF 


a  TURF  Naturally  3D  Cartesian 
a  Considered  2D-Axisymmetric  TURF 

a  Simple  DSMC  but  Merge  Complex 
(Conservation  on  V||Vj_ -Quadtrees?) 


Standard  3D  DSMC 


Robert  Martin  (AFRL/RQRS)  distribution 
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o  TURF  Naturally  3D  Cartesian 
o  Considered  2D-Axisymmetric  TURF 

o  Simple  DSMC  but  Merge  Complex 
(Conservation  on  V||V^-Quadtrees?) 

o  Opted  to  Run  Coarse  Full  3D 
(Simplified  Boundary  Conditions) 

o  3D  Expensive  at  Tractable  Resolution 
o  Added  Collision  Sub-Cells  to  TURF 
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Challenges  for  Jet  in  TURF 


a  TURF  Naturally  3D  Cartesian 
a  Considered  2D-Axisymmetric  TURF 

a  Simple  DSMC  but  Merge  Complex 
(Conservation  on  V||Vj_ -Quadtrees?) 

a  Opted  to  Run  Coarse  Full  3D 
(Simplified  Boundary  Conditions) 

a  3D  Expensive  at  Tractable  Resolution 
a  Added  Collision  Sub-Cells  to  TURF 
a  Fractional  DSMC  Controls  Np/Cell 
a  Linear  Density  Obscures  Results 
a  Issue  Clearer  with  NpAr/Cell 


Standard  3D  DSMC 


Multi-Species  Fractional 
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Challenges  for  Jet  in  TURF 


a  TURF  Naturally  3D  Cartesian 
a  Considered  2D-Axisymmetric  TURF 

a  Simple  DSMC  but  Merge  Complex 
(Conservation  on  V||Vj_ -Quadtrees?) 

a  Opted  to  Run  Coarse  Full  3D 
(Simplified  Boundary  Conditions) 

a  3D  Expensive  at  Tractable  Resolution 
a  Added  Collision  Sub-Cells  to  TURF 
a  Fractional  DSMC  Controls  Np/Cell 
a  Linear  Density  Obscures  Results 
a  Issue  Clearer  with  NpAr/Cell 

XAr  to  RX-Plane  — >  for  Detailed  Results... 
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Preliminary  TURF 


9  Standard  DSMC  Poor  Results 
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Preliminary  TURF  Results:  He:Ar  Jet 
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O  Standard  DSMC  Poor  Results 
Q  2x2x2  Collision  Cell  Improves  Standard  DSMC 
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Preliminary  TURF  Results:  He:Ar  Jet 


Q  Standard  DSMC  Poor  Results 
Q  2x2x2  Collision  Cell  Improves  Standard  DSMC 
Q  SWPM+Octree  Significantly  Better  (2x2x2  Collision  Cell) 


Future  Directions: 


The  df- Boltzmann  for /  —  feq  -\-fdev 

o  df-  Concept  Old 

o  How  to  make  df  cheaper  than  full-/ 
(Must  adapt  DOF  Usage..?) 
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Future  Directions: 


The  df- Boltzmann  for /  —  feq  -\-fdev 

o  df-  Concept  Old 

o  How  to  make  df  cheaper  than  full-/ 

(Must  adapt  DOF  Usage..?) 

o  Recent  Progress  using  ±£-weight  Particles 
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Future  Directions: 
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o  Recent  Progress  using  ±£-weight  Particles 
9  Requires  Remapping  due  to  Particle  Growth 
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Future  Directions: 
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o  df  -  Concept  Old 

o  How  to  make  df  cheaper  than  full-/ 

(Must  adapt  DOF  Usage..?) 

o  Recent  Progress  using  =b£-weight  Particles 
9  Requires  Remapping  due  to  Particle  Growth 

o  SWPM  has  Similar  Issue 

(Basis  for  Octree  B2B  Collisions) 
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The  df- Boltzmann  for /  —  feq  -\-fdev 

o  df-  Concept  Old 

o  How  to  make  df  cheaper  than  full-/ 

(Must  adapt  DOF  Usage..?) 

o  Recent  Progress  using  ±£-weight  Particles 
9  Requires  Remapping  due  to  Particle  Growth 

®  SWPM  has  Similar  Issue 

(Basis  for  Octree  B2B  Collisions) 

o  Octree:  Bin  Moments,  Mi  — >  Particles 
o  Root  Bin  Sum  — y  Equilibrium  Moments 


Robert  Martin  (AFRL/RQRS) 


DlSTRIBU] 


Octree  df- Boltzmann 
Bin-to-Bin  Collisions 


_ / 

'O  <\  O' 


5E;  DlST 


19/21 


Future  Directions: 


The  df- Boltzmann  for /  —  feq  -\-fdev 

o  df-  Concept  Old 

o  How  to  make  df  cheaper  than  full-/ 
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Sample  Collisions  using 
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Current  Results: 


©  SWPM+Octree  Option  for  Variable  Weight  Mixture  Collisions 
o  Multiple  Octree  Merge  only  Modification  for  Multi-Species 
©  Initial  Verification  vs.  Standard  Shock  Cases 
©  Merge/Target  Enables  Direct  Control  of  Noise 
©  Unidentified  Systematic  Error  with  1 : 1  Target 
©  Initial  Testing  on  3D  Mixture  Expansion  Better  with  SWPM+Octree 

Future  Efforts: 

©  Additional  Investigation  of  Error  Source  for  Disparate  Weights 
©  Improved  Merge/Control  of  Tails 
©  Apply  to  Reacting  Flow 
©  Adaptation  for  5f 


Thank  You 


This  material  is  based  upon  work  supported  by  the  Air  Force  Office  of 
Scientific  Research  under  award  number  FA9550-17RQCOR465. 

Any  opinions,  finding,  and  conclusion  or  recommendations  expressed  in  this  material  are  those  of  the  author  and 
do  not  necessarily  reflect  the  views  of  the  United  States  Air  Force. 
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